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lO ! Abstract 

We analyze the statistical mechanical properties of n-detectors in arbitrary states 
O^' of motion interacting with each other via a quantum field. We use the open system 

^ ! concept and the influence functional method to calculate the influence of quantum fields 

on detectors in motion, and the mutual influence of detectors via fields. We discuss 
the difference between self and mutual impedance, advanced and retarded noise, and 
^ , the relation of noise-correlations and dissipation-propagation. The mutual effects of 

5^ \ detectors on each other can be studied from the Langevin equations derived from 

the influence functional, as it contains the backreaction of the field on the system 
self-consistently. We show the existence of general fluctuation-dissipation relations, 
and for trajectories without event horizons, correlation-propagation relations, which 
succinctly encapsulate these quantum statistical phenomena. These findings serve to 
clarify some existing confusions on the accelerated detector problem. The general 
methodology presented here could also serve as a platform to explore the quantum 
statistical properties of particles and fields, with practical applications in atomic and 
optical physics problems. 
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1 Introduction 



The physics of accelerated detectors became an interesting subject of investigation when 
Unruh showed that a detector moving with uniform acceleration sees the vacuum state of 
some quantum field in Minkowski space as a thermal bath with temperature Tu = ha/2TTckB 
This seminal work which uses the structure of quantum field theory of Rindler space 
explored by Fulling 0, and the Bogolubov transformation ideas invented by Parker for cos- 
mological particle creation ^ earlier, draws a clear parallel with the fundamental discovery 
of Hawking radiation in black holes The discovery of Unruh effect (see also Davies ^) 
sets off the first wave of activities on this subject. The state-of-the-art understanding of the 
physics of this problem in this first stage of work is represented by the paper of Unruh and 
Wald 0. We refer the readers to the reviews of Sciama, Candelas and Deutsch 0, Tagaki 
g] and Ginzburg and Frolov 

The second stage of investigation on this problem was initiated by the inquiry of Grove 
10| , who challenged the prevailing view and asked the question whether the detector actually 



radiates. This was answered in the negative by an inspiring paper of Raine, Sciama and Grove 
||TT| (RSG henceforth) who considered an exactly solvable harmonic oscillator detector model 
and analyzed what an inertial observer sees in the forward light cone of the accelerating 
detector via a Langevin equation. Unruh [1^ performed an independent calculation and 
concurred with the findings of RSG to the extent that the energy-momentum tensor of the 
field as modified by the presence of the accelerating detector vanishes over most of the 
spacetime (except on the horizons). However, he also showed the existence of extra terms 
in the two-point function of the field beyond its value in the absence of the accelerating 
detector, and argued that these terms would contribute to the excitation of a detector placed 
in the forward light cone. These terms were missed out in RSG. Following these exchanges, 
there was a recent renewed interest in this problem, notably the series of papers by Massar, 
Parantani and Brout (MPB), who gave a detailed analysis via Hamiltonian quantum 
mechanics of the two-point function and pointed out that the missing terms contribute to 
a polarization cloud around the accelerating detector; Hinterleitner [Q, who independently 
discussed the backreaction of the detector on the field using a slightly different yet exactly 
solvable model and arrived at similar conclusions to MPB; and Audretsch and Miiller ||15|| , 
who explored nonlocal pair correlations in accelerated detectors. However, the physical 
significance of the polarization cloud, its connection to the noise experienced by another 
detector, and to the inherent correlations in the free Minkowski vacuum, remain largely 
unexplored. 

Beginning with this work we would like to add a new dimension to this problem and 
open up the third stage of investigation. The new emphasis is in exploring the statistical 
mechanics of particles and fields, and in particular, moving detectors on arbitrary trajec- 
tories. We analyze the stochastic properties of quantum fields and discuss this problem in 
terms of quantum noise, correlation and dissipation. We use the open system concept and 
the influence functional method to treat a system of n detectors interacting with a scalar 
field. This method enables one to examine the influence of detectors in motion on quantum 
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fields, the mutual influences of detectors via flelds, as well as the backreaction of flelds on 
detectors in a self-consistent manner. 

As explained earlier [1^, the influence functional method is a generalization of the power- 
ful effective action method in quantum fleld theory for treating backreaction problems, which 
also incorporates statistical mechanics notions such as noise, fluctuations, decoherence and 



dissipation. Indeed, one of us has long held the viewpoint that IT^, |T8[) to get a more 
profound understanding of the meaning of Unruh and Hawking effects and the black hole 
information and backreaction problems one cannot be satisfled with the equilibrium ther- 
modynamics description. It is necessary to probe deeper into the statistical properties of 
quantum flelds, their correlations and dynamics, coherence and decoherence of the particle- 
fleld system, the relation of quantum noise and thermal radiance, fluctuation-dissipation 
relations, etc. Earlier investigation of correlation and dissipation in the Boltzmann-BBGKY 
scheme [19, 21] and the properties of noise and fluctuations in the Langevin framework 
P^ , are essential preparations for tackling such problems at a deeper level. [] 

With this theoretical perspective in mind, we have recently begun a systematic study 
of the accelerated detector problem p5| , We show that thermal radiance can be un- 
derstood as originating from quantum noise under different kinematical (moving detector) 
and dynamical (cosmology) excitations. The aim of this paper is to 1) show on both the 
conceptual and technical levels the power and versatility of this new method, 2) settle some 
open questions and clarify some existing confusions, such as the existence of radiation and 
polarization, solely from an analysis of detector response, 3) introduce new concepts such 
as self and mutual impedance, advanced and retarded noise, fluctuation-dissipation and 
correlation-propagation relations using the accelerated detectors problem as example, and 
flnally 4) suggest new avenues of investigations into the statistical mechanics of particle and 
flelds, including black hole physics. 

Employing a set of coupled stochastic equations for the detector dynamics, we analyze 
the influence of an accelerated detector on a probe which is not allowed to causally influence 
the accelerated detector itself. We flnd, as did |jTl], |12|, ^ that most of the terms in the 
correlations of the stochastic force acting on the probe cancel each other. This cancellation is 
understood in the light of a correlation-propagation relation, derived as a simple construction 
from the fluctuation-dissipation relation for the accelerated detector. Such a relation can 
be equivalently viewed as a construction of the free fleld two-point function for each point 
on either trajectory from the two-point function along the uniformly accelerated trajectory 
alone. The remaining terms, which contribute to the excitation of the probe, are shown to 
represent correlations of the free fleld across the future horizon of the accelerating detector. In 
this problem, the dissipative properties of either detector remain unchanged by the presence 
of the other. This happens because the probe cannot influence the accelerated detector. 
However, the stochastic force acting on the probe plays a non-trivial role. 

We also consider the problem of two inertial detectors which can backreact on each 
other. This mutual backreaction changes the self-impedance functions of these detectors, 
and introduces mutual impedances as well. The dissipative properties of each detector 

^Our view is thus most akin to that espoused by Sciama p^. 
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are thus altered due to the presence of the other one. This physical effect is in a sense 
complementary to the effects manifested in the accelerated detector problem, where the 
probe does not backreact on the accelerated detector. 

The paper is organized as follows. In Section 2 we develop the influence functional 
formalism describing the influence of a massless scalar fleld on a system of an arbitrary 
number of detectors moving on arbitrary trajectories. The fleld modes are integrated out 
in this formalism, and effective stochastic equations of motion for the various detectors are 
obtained. In Section 3 we consider some applications of this formalism to three simple cases, 
the primary one being the analysis of two inertial detectors coupled to the same quantum 
fleld. In Section 4 we treat the RSG excitation of a probe in the presence of a uniformly 
accelerating detector. In Section 5 we show the existence of fluctuation-dissipation relations 
governing the detector system. These relations are used as a starting point for obtaining more 
general relations between the correlations of various detectors and the radiation mediated 
by them. Such relations are also discussed in the speciflc context of the RSG model. Finally, 
in the Appendix, we point out problems associated with the uncorrelated detector-fleld 
initial state in a minimally coupled model, and argue that these problems are removed in a 
derivative coupling model. We present a simple prescription for switching from one model 
to the other. 

2 Scalar Electrodynamics or Minimal Coupling Model 

The paper by Raine, Sciama and Grove uses the scalar electrodynamic or "minimal" coupling 
of oscillators to a scalar fleld in 1+1 dimensions. This couphng provides a positive deflnite 
Hamiltonian, and is of interest because it resembles the actual coupling of charged particles 
to an electromagnetic field. In this section, we derive the infiuence functional describing 
the effect of a scalar field on the dynamics of an arbitrary number of detectors modelled as 
minimally coupled oscillators. The detectors move along arbitrary trajectories. We assume 
that the fleld and the system of detectors are initially decoupled from each other, and that 
the field is initially in the Minkowski vacuum state. The formalism can be simply extended 
to higher dimensions, and to different choices of initial state for the fleld. We also obtain 
coupled Langevin equations for the detector system. 

2.1 Influence functional for N arbitrarily moving detectors 

Consider detectors i = l,..iVinl + l dimensions with internal oscillator coordinates Qi{Ti), 
and trajectories {xi{Ti),ti{Ti)) , Tj being a parameter along the trajectory of detector i. In 
the following analysis, we do not need to assume that Ti is the proper time, although this is, 
in most cases, a convenient choice. However, we will assume hereafter that the trajectories 
{ti{Ti),Xi{Ti)) are smooth and that the parameters are chosen such that ti{Ti) is a strictly 
increasing function of Tj. 
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The detectors are coupled to a massless scalar field t) via the interaction action 



Y,^if' driSi{ri)—^(f){xi{ri),ti{ri)). (2.1) 



Here, T is a global Minkowski time coordinate which defines a spacelike hypersurface, Cj 
denotes the coupling constant of detector i to the field, Sj(rj) is the switching function for 
detector i (typically a step function), and is the inverse function of ti. t1^{T) is therefore 
the value of at the point of intersection of the spacelike hypersurface defined by T with 
the trajectory of detector i. Note that the strictly increasing property of ti{Ti) implies that 
the inverse, if it exists, is unique. 

The action of the system of detectors is 



"^osc 2 



1 

r dn[{dr,Qif-^m- (2-2) 

The scalar field action is given by 

S field = \ dt I dx[{dt<Pf - {d,<Pf] (2.3) 

and the complete action 

S — S field + Sosc + Sint- (2.4) 

Expanding the field in normal modes, 

(f){x, t) = Jy JZi'ikit) cos kx + q^{t) sin kx] (2.5) 

' k 

where Yl'k denotes that the summation is restricted to the upper half k space, k > 0. Then 
the action for the scalar field is given by (cr = +, —) 

Sf^eld = ltMf-u^lql] (2.6) 

k,(7 

and the interaction action is 

Sint = Z^^iW- / dTi——X 



oo dTi 

Yl\.^k(ti(n)) COSkXi{Ti) + q^{ti{Ti)) SVa.kXi{Ti)\Si{Ti) 
k 



poo nt^ {1 ) 

- dt I dTi5(t - ti{Ti)) 

ij J —CO J —oo 



dTi 



X 



Y,[q]^(t) COS kxi(Ti) +q^(t) sin kxi(Ti)]si(ri). (2.7) 

k 
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We have tj(rj) < T, which follows from Tj < t^^(T) and the property that U^Ti) is a strictly 
increasing function. Hence we may replace the upper limit of the dt integration by T. This 
replacement leads to the expression: 

S^nt = ~i: r dtJl{t)ql{t) (2.8) 

^ j^^ dn5{t - tM))-^ul{n)s,{n) (2.9) 



where 



and 

ul{T.i) = coskxi^Ti); u^i^i) = sinkxi^Ti). (2.10) 

The action Sfieid+Smt therefore describes a system of decoupled harmonic oscillators each 
driven by separate source terms. The zero temperature influence functional (corresponding 
to the initial state of the field being the Minkowski vacuum) for this system has the form 



k,cr 

where 



^[J,J']=exp{-i^ r ds r dsyas)-j',^mCkis,s')r,is')-Cds,sV'k''is')]} (2.11) 



C, = u, + tf^, = ^e-'^'=^'-''\ (2.12) 

If the field is initially in a thermal state, the influence functional has the same form as above, 
and the quantity (k becomes 

Ck = cosa;fc(s - s') -ismujk{s - s')], (2.13) 

f5 being the inverse temperature. We shall restrict our attention to the zero temperature 
case. 

Substituting for the J|!"'s in the influence functional, and carrying out the 5-function 
integrations, one obtains 



^[{QrAQ}] = exp--{}^y_ dnsiir,) dTjSj{Tj)[— - —\ X 

where 

Zi,{ri,r;) = |e,e,^a(t,(r,),t,(r;))<(r,)<(rj). (2.15) 
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In the above, the continuum hmit in the mode sum is recovered through the replacement 
— > ^ Jo° dk. We then obtain, after substituting for and (k, 

Z^An.r'^) = ^ |^°°^e-^'=(*'(-)-*^(-i»cosMx.(r,) -x,(r;)). (2.16) 

In this form, Zij is proportional to the two point function of the free scalar field in the 
Minkowski vacuum, evaluated for the two points lying on trajectories i and j of the detector 
system. It obeys the symmetry relation 

Z,,(r„r;)=Z;(r;,r,) (2.17) 

Corresponding to (2.12), we may also split Zij into its real and imaginary parts. Thus 
we define 

Zijin, Tj) = Uijin, Tj) + ip,ij{Ti, Tj) (2.18) 



where 



J^ij{Ti, Tj) = ^ ^ ^ COS k{ti{Ti) - tj{Tj)) COS k{Xi{Ti) - Xj{Tj)) 



g -g . POO fj^J^ 

fiij{Ti,Tj) = — ^ — sin k{ti{Ti) -tj{Tj)) cos k{xi{Ti) - Xj{Tj)). (2.19) 

i) and fl are proportional to the anticommutator and the commutator of the field in the 
Minkowski vacuum, respectively. 

The quantities Z^j are also conveniently expressed in terms of advanced and retarded null 
coordinates Viin) = U{Ti) + Xiln) and Uiln) = Uiri) - Xi{Ti), as 

z,,in,T^) = z^in,T^) + zrin,T^) (2.20) 

where 

Z^^in,T') = ^ r ^^e-*^("'(^^)-"^(^i« (2.21) 
An Jo k 

and the superscripts a and r denote advanced and retarded respectively.^ Similar decompo- 
sitions for Uij and flij thus follow. 

The influence functional, along with the action for the detector system, can be employed 
to obtain the propogator for the density matrix of the system of detectors. This propogator 
will contain complete information about the dynamics of the detectors. However, we shall 
take the alternative approach of deriving Langevin equations for the detector system in order 
to describe its dynamics. 

^ The terminology 'advanced' and 'retarded' refers to the nuU coordinates. Equivalently, they can be 
called 'left-moving' and 'right-moving', respectively, when the sense of motion refers to the future direction 
in time. This terminology is used in wave theory and string theory. 
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2.2 Langevin equations 

In this subsection, we wish to derive the effective stochastic equations of motion for the N- 
detector system. In the previous subsection, we integrated out the field degrees of freedom. 
The effect of this is to introduce long-range interactions between the various detectors. 

Going back to the form (2.11) for the influence functional, we define the centre of mass 
and relative variables 

j^is) = ms)+Ji^{s))/2 

Jk^is) = r,{s)-Jl^{s). (2.22) 
Correspondingly, we also find it convenient to define 

Qtin) = mn) + Q[{n))/2 

gr(r,) = Q.(r,)-Q[{n). (2.23) 



Then Equation (2.11) yields 

lA 



^[J,J']| = eM-TT.1 ds'J^'^{s)vk{s,s')J^'^{s')} (2.24) 

- / KAWuPru]) exp ~ E dsJ;'^{s)il{s). (2.25) 



I ^ I is the absolute value of ^, containing the kernel v^. The phase of T contains the kernel 
Hk. In the second equality, we have used a functional gaussian integral identity, P[^%] being 
the positive definite measure 

Pm = iVexp{-^ r ds r ds'Ck{sH\s,m{s')} (2.26) 

normalized to unity. It can therefore be interpreted as a probability distribution over the 
function space 

The influence functional can thus be expressed as 

n{Q},{Q'}] = <exp{~j2 r dsJr{sm{s)+2 r ds'f^k{s,s')Jr{s')]}> 

= <exp^Sinf> (2.27) 

where < > denotes expectation value with respect to the joint distribution n'^^P[^^]. Sinf 
will be called the stochastic influence action. We find 

<Ckis)> = 0, 
<{^ki^),^k'i^')}> = h6kk'Saa"^k{s,s') (2.28) 
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where { , } denotes the anticommutator. 

Substituting for J'j^'^ and J^°^ in terms of the detector degrees of freedom {Qj}, the 
stochastic influence action Smj is obtained as 



S^nf = -Y.^ dn-^s,{n)[r^,{n)+2Y.j' d^-^s.ij'^f^^.ir^y^ (2.29) 

with 

V^{r^) = e, ^ <(r.)4"(t.(r,)). (2.30) 

From Equation (2.29) we see that the quantities p^ij, i j mediate long-range interactions 
between the various detectors and the quantities flu describe self-interaction of each detector 
due to its interaction with the field. This self-interaction typically manifests itself as a 
dissipative (or radiation reaction) force in the dynamics of the detectors. We will, therefore, 
refer to /2jj, i ^ j as a "propagation kernel", and "dissipation kernel" . 

We now turn to the interpretation of the quantities rji. They appear as source terms in 
the effective action of the detector system. Also, being linear combinations of the quantities 

they are stochastic in nature. Indeed, from Equations (2.28) and (2.30) we can obtain 



< Vii^t) > = 0, 

< Mr,),v,{r;)} > = e,e,Y.jz ul{n)ut,{T'^){j) < 4^(t.(r,))4<(t,(rj)) > 



= /i%(r,,r;). (2.31) 

Thus % appears as a correlator of the stochastic forces rji and rjj. Along a fixed trajectory, 
this correlation manifests as noise in the detector dynamics. Hence we call vu a "noise 
kernel" and z/jj, i ^ a "correlation kernel" .0 

The full stochastic effective action for the A^-detector system is given by 

•^e// = Sosc + Sinj (2.32) 

We may now express this in terms of the variables Ql and defined earlier. Thus we 
obtain 

Seff = I dT,[Q;Qt -n^Q^Qt -Q7s,{T,Mn) 

-2q-s,{t,)J2 i ^ ^ dT;Qp,{T;)fi,,{T,,r;)] (2.33) 



•^The distinction between noise and correlation is unnecessary from the point of view of the field. 'Noise', 
as used here, also represents free field correlations for points on a single trajectory. However, from the point 
of view of each detector, these two quantities play a different role. Hence the choice of terminology. 
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where /^ ^ f , ^ ^. 

Extremizing the effective action with respect to Ql and setting Qi = Q'^ at the end |p2|| , 
we obtain a set of coupled equations of motion, the Langevin equations, for the system of 
detectors: 

^ - 2Ey_'^'^*'^""cir;.,<rj)A(,^(^^)^^^.(^^,^;))^ + ^^2^ ^ ^(..(r,)r^.(r,)). (2.34) 

Due to the back-reaction of each detector on the field, and consequently on other detectors, 
the effective dynamics of the detector system is highly non-trivial and, as such, can be solved 
in closed form only for simple trajectories or under simplifying assumptions such as ignoring 
the back-reaction of certain detectors on the field. For instance, if we choose to ignore the 
back-reaction of detector i on the field, this can be effected by setting jlji = 0, for all j, 
including j = z, while at the same time keeping fiij ^ for j ^ i. The particular case 
jlii = amounts to ignoring the radiation reaction of detector i. This is necessary because 
the radiation reaction effect arises due to a modification of the field in the vicinity of the 
detector as a consequence of the back-reaction of the detector on the field. 

Of course, it is in general inconsistent to ignore the back-reaction of a detector, as it leads 
to a direct violation of the symmetry (2.17). As is well-known, it also leads to unphysical 
predictions. For example, in the treatment of an atom on an inertial trajectory, coupled to a 
quantum field, balance of vacuum fluctuations and radiation reaction is necessary to ensure 
the stability of the ground state. As explained above, ignoring back-reaction implies ignoring 
the radiation reaction force. Such a treatment would render the ground state unstable. 

However, in certain cases, the quantities ftji may not contribute to the dynamics of 
detector j, as in Section 4 below, where the trajectory of one detector is always outside the 
causal future of the other one. Hence there is no retarded effect of one of the detectors on 
the other. 

Our formal treatment of the detector-field system is exact in that it includes the full back- 
reaction of the detectors on the field, which is manifested in the coupled Langevin equations 
of the various detectors. The coupled equations of motion give rise to a sort of "dynamical 
correlation" between the various detectors. Non-dynamical correlations also occur because 
of the intrinsic correlations in the state of the field (Minkowski vacuum). These correlations 
are purely quantum-mechanical in origin, and they are reflected in the correlators of the 
stochastic forces, Uij. Correlations between stochastic forces on different detectors induce 
correlations between the coordinates Qi of different detectors. 

As we shall show in a later section, our exact treatment makes it possible to demonstrate 
the existence of generalized fluctuation-dissipation and correlation-propagation relations gov- 
erning the detector system. 
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3 Examples 



In this and the following section, we consider some applications of the Langevin equations 
derived in the previous section to the cases of a single detector in the Minkowski vacuum 
moving on an inertial trajectory, a single detector on a uniformly accelerated trajectory, two 
detectors on inertial trajectories, and the case of one detector on a uniformly accelerated 
trajectory and another one on an arbitrary trajectory, functioning as a probe. 

The first two examples serve to illustrate the formalism, and describe the well-known 
physical effects of the dressing of a particle by the field and the thermal Unruh noise expe- 
rienced by a uniformly accelerated particle. 

In the example of two inertial detectors, we introduce the notions of "self" and "mutual" 
impedance which govern the response of either detector. The effect of the back-reaction of 
each detector on the field and consequently on the other detector is to introduce the so-called 
mutual impedance in the detector response as well as to modify the self-impedance of each 
detector from its value in the absence of the other one. 

In the next section we shall consider the example of one detector on a uniformly acceler- 
ated trajectory and a probe, which moves along an unspecified trajectory. Wc switch on the 
probe after it intersects the future horizon of the uniformly accelerated detector, so that it 
cannot causally influence the uniformly accelerated one. Thus the uniformly accelerated de- 
tector in this case is effectively in an unperturbed Unruh heat bath, and this situation mimics 
most closely the RSG model. The missing terms in the RSG analysis, which contribute to 
a polarization cloud around the accelerated oscillator, but not to the energy momentum 
tensor, lead to a modified noise kernel in the Langevin equation for the probe. 

In all cases, we can solve exactly for the detector coordinates, at least in the late time 
limit (this limit is actually realized at any finite time t ^ — oo when the two detectors have 
been switched on forever, and corresponds to the neglect of transients in the solutions for 
the detector coordinates). 

3.1 One inertial detector 

Consider the case of one detector moving on an inertial trajectory x(r) = 0, i(r) = r, and 
switched on forever (s(r) — 1). The noise and dissipation kernels take the form 




(3.2) 



(3.1) 



The Langevin equation becomes 



d^Q dQ 



drj 
d^ 



(3.3) 



with 




(3.4) 



11 



It will be convenient to define the dissipation constant ^ — \ ■ We will restrict our attention 

to the undcrdamped case (7 < f^o )■ 
Introducing the Fourier transform 

/oo 
-00 

Q{uj) ^ ^ r due-^'^Qir) (3.5) 
and similarly for ?7(t), we obtain 

Qiuj) = Xujfi{uj) (3.6) 
with the impedance function Xw defined as 

Xu^^iuj{-uj'^ + ^l + 2iuj^)-\ (3.7) 

In the above solution for the detector coordinate in frequency space, it should be noted 
that transients have already been neglected. Transient terms correspond to delta functions 
in frequency space, the coefficients of these delta functions being determined by the initial 
conditions. For the complete solution these terms should be added to the right hand side of 
Equation (3.6). We may thus obtain 

< {Q{uj),Q{uj')} >= x-Xo;' < {fi{uj)Muj')} > (3.8) 

where 

<{fj{u;),fj{u;')}> = —j^Jr j_jT'e-^-^e-^^ <{r^{T),r^{T')}> 

= h^5{uj + uj')[e{uj) - e{-uj)\. (3.9) 

TIOJ 

We can therefore obtain the correlator of Q{t) and Q{t'), as 

/oo fOO . 

duj / du'e^'e'^ ' X.X.' < > 
-00 J—oo 

2^7 dk 



r^lXfeTcosMr-r'). (3.10) 
JO k 



71 Jo k 

3.2 One accelerated detector: Unruh effect 

In the case of an accelerated detector moving on the trajectory x{t) — cosh ar, t(r) = 
sinh or, and s(r) = 1 (r being the proper time along the accelerated trajectory), the 
noise and dissipation kernels take the form: 



v{t,t') = — [ — COS —(sinh ar — sinh ar') COS —(cosh ar — cosh ar') (3.11) 
2tt Jo k a a 

il(T,r') — — — J — sin —(sinh or — sinh ar') COS -(cosh ar — cosh ar'). (3.12) 
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These kernels can be decomposed into advanced and retarded parts, by writing, for example 



/■°° dk 



y'^ + y^ = I 



COS -(e"" - e""') + cos -(e"'^" - e 
a a 



(3.13) 



An Jo k 

We can then use the changes of variables k — > |e^2*-'^~'^''' to obtain 

showing that the noise felt by the accelerating detector is isotropic. 

One can also make a similar simplification for the kernel fi. These expressions can then 
be further simplified p5| , p3| by means of the integral transform p7 



2 roo 

ik sinh h 



IT Jo 

and the formula 



da cos{ba + iTTa/2)Kia{k) (3.15) 



where Kia{a) is a Bessel function of imaginary argument, to yield 

z>(r, r') = — / — coth( — )cosA;(r — r') (3-17) 
ZTT Jo k a 

Mr,r') = --y^ ^sinA:(r-r'). (3.18) 

The noise experienced by the detector is thus stationary and the factor coth(— ) in the noise 
kernel shows that it is also thermal, at the Unruh temperature fc^T = |^ (we have chosen 
units such that c = 1). 

The dissipation kernel remains identical to that of the inertial detector. Based on the 
property that the two-point function of a free field on an accelerated trajectory evaluated in 
the Minkowski vacuum state is identical to the two-point function on an inertial trajectory 
evaluated in a thermal state at the Unruh temperature, this fact can be explained as follows: 
The dissipation kernel is proportional to the commutator of the free quantum field evaluated 
in whatever state the field is in. However, the commutator of a free field for any two points 
is just a c-number, hence its expectation value is independent of the state of the field. In 
particular, it does not distinguish between a zero temperature and a thermal state. So the 
dissipation kernel is identical to that in the inertial case. The anticommutator is, however, 
an operator whose expectation value depends on the state of the field, and therefore shows 
the familiar departure from the inertial case. 

The Langevin equation for the detector coordinate is 

d^Q e'^dQ , df] 
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with 

<Mr),v{r')}>^hiy{T,T'). (3.20) 
Similar to the inertial detector case, we find 



Qiuj)=x^viuj). (3.21) 



Also 



1 POD rOD . 



= ;i^coth(— )% + u;')- (3-22) 

TTCU a 

Combining the two equations, 

/oo roo . 
do; / du'e^'e'^ ' x.X.' < {fl{u;),fj{u;')} > 

-oo J —oo 

= — - -r Xk coth — cos/c r - r' 3.23) 

n Jo k a 

with the impedance function Xk as defined in the inertial case. 

3.3 Two inertial detectors: self and mutual impedance 

We now consider the case of two detectors moving on the inertial trajectories Xi{ti) = —Xo/2, 
X2{t2) = xq/2 and ti(ri) = t2(T2) = t, coupled to a scalar field initially in the Minkowski 
vacuum state, with coupling constants ei,2. They are separated by a fixed coordinate distance 
xq. As before, we will assume that both detectors have been forever switched on, i.e. ^^(t) = 

1,2 = 1,2. 

It will be convenient to express the noise, dissipation, correlation and propagation kernels 
as the real and imaginary parts of the functions defined earlier. Then, for the two-detector 
system, we obtain 



Zi2(T,T') = i?2i(T,T') = ^jf°°^e-«<'--''cosfa„. (3.26) 
The coupled Langevin equations for the system are 

d^Qi , efdQi , 6162 rfgs I .^2n o7^ 

^+2"^ + ^^'---+^^^^ - d7 ^^-^^^ 

d^Q2 el dQ2 , 6162 dQi ^ di]2 , . 
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where t — xq is the retarded time between the two trajectories, and 

< Mr),Vj{r')} >= hu,,{T - t'). (3.29) 

As before, we define 71 2 = and introduce Fourier transforms to obtain the corre- 
sponding equations in frequency space. Then we obtain 

QiH = xL'^^i(^) - xi,^^Q2(a;) (3.30) 

^2(0;) = xi?)r/2(a;)-^e--°xL'^Qi(a;) (3.31) 

where 

xL'^'^'^ = icui-u^ + n?,2 + 2iiu-fi^2y\ (3.32) 

The functions xL^'*'*'^'' of course, what the impedance of each detector would be in the 
absence of the other one. However, the effect of introducing a second detector is, as we 
shall see, to modify the "self - impedance" of each detector as well as introduce a "mutual 
impedance" which describes, for instance, the response of detector 1 to the force 772- 
Indeed, plugging the equation for Qi in the equation for Q2, we have 

Q2{uj) = L22{uj)fj2{uj) + L2l{Lu)fji{uj) (3.33) 

where L22 is the modified self - impedance of detector 2 due to the presence of detector 1, 
and L21 is the mutual impedance: 

L22(^) = xL'Hl-47i72e-^^-^°xL')xL'))-' 

L2iiu;) = -2y^e— °xL'^xL'ni-47i72e-^^--°xL'^xL'^)-^ (3.34) 

The impedances Lu and L12 and the corresponding equation for Qi are obtained by an 
interchange of indices 1 and 2 in the above equations. 
We note the symmetry 

L21 — Li2- (3.35) 

The correlator < {Qi{ijj),Qj{ijj')} >, i,j — 1,2 is therefore obtained from equation (3.32) 
and its counterpart, as 

2 2 

< {Q,{uj),Qj{uj')} >= Y^Y.^i-H^A'^') < fe(^),^/3(^')} > • (3.36) 

a=l 0=1 

The above equation is to be viewed as a generalization of (3.8) to the two-detector case. 

Suppose we now wish to solve for the correlator of Q2- Then, taking Fourier transforms 
as before and simplifying, 

duj / dJe'^^e'^^ < {Q2{u;),Q2{u;')} > (3.37) 
-00 J —00 

2/172 f°° dk 2 n , fA I (1) |2 (1) 

= y I ^22(A:) I [1 + 71(471 \Xk \ -Xk -Xk ) 

-7i(xi'^e-^^'^^o+^(i)*e2.^.o)]. 
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The second term in the square brackets vanishes as a consequence of the identity Xk^+X^k^* — 
^71 I X^k^ which is a form of the fluctuation dissipation relation for detector 1. The 
remaining terms simplify to yield 

<{Q2{t),Q2{t')}>^ (3.38) 

^ ? I L22{k) r [1 - 271 I x« r (271 cos2fca;o + ^^^^ sin 2fca;o)] x 
COS k{T — t'). 

As before, the correlator of Qi is obtained by interchanging the indices 1 and 2 in the 
above equation. 



4 Uniformly accelerated detector and probe: advanced 
and retarded noises 

In this section, we first consider two detectors on arbitrary trajectories, with the following 
constraints: a) both trajectories are everywhere timelike, b) one of the trajectories possesses 
past and future event horizons, which are chosen to be the null lines v — Q and u — Q 
respectively, c) the detector on the other trajectory is switched on at u = and this trajectory 
does not possess a future horizon. 

Because of constraint c) , the second detector cannot causally influence the first one, and 
thus it functions as a probe in the field modified by the first detector. 

Later in the analysis, we shall specify the trajectory of the detector with horizons as being 
a uniformly accelerated one. We shall continue to assume that the probe cannot causally 
influence the uniformly accelerated detector by means of the switching condition. If it were 
allowed to do so, this would lead to a deviation of the noise experienced by the uniformly 
accelerated detector from the precise thermal form. 

We will label the detector with horizons as detector 1 and the probe as detector 2. The 
switching condition S2(t2) = 0{u2{t2)) for the probe leads to a closed Langevin equation for 
detector 1: 

(fQi el dQi _ u i ^ 

^+ 2"d^ + ^^^^-d^- ^^-^^ 
This is just a consequence of the fact that the trajectory of detector 1 hes outside the causal 
future of the probe. The arguments which lead to the above local form of dissipation or 

radiation reaction for a general timelike trajectory are outlined in the next section. 

Introducing Fourier transforms and the impedance functions ''''^^^ ^ defined earlier, we 
have, 

Q,iu) = x^^^,{u). (4.2) 
Consider now the Langevin equation for detector 2: 



dT' 



2 



' / C?T-2 — (e(li2(T2))/i22(T2,T^))-rT 
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J-oo UT2 CLTi aT2 

We find 



^{~^22{t2A)) = -272<^(t2-t^) (4.4) 

^(/i2i(r2,rO) = -V7l7^[^5(^2(T2)-^i(rO) + ^5(«2(r2)-«i(rO)] (4.5) 

where 71,2 are defined as in the two inertial detector case. The second term in (4.5) vanishes 
identically because U2{t2) > and ui{ti) < (li = is a future horizon for detector 1). 
Since = is a past horizon for detector 1, we have vi{t[) > and the first term simplifies 
to yield 

-^{il2l{ru t',)) = -^2^5{tr - T[)e{v2{T2)) (4.6) 
dT2 OT2 

where we have defined the retarded time tr = f r^(^2(T2)). This is well-defined since it occurs 
only in expressions in which i'2(t2) > 0. 

Thus we obtain the dynamical equation for the probe, which depends, as expected, on 

Qi: 

(fQ2 , e| dQ2 , „2^ 6162 (Itr dQi / x ^ n ^\ 



Consider the quantity 

dr]2 6162 dTR dQi 

which is a source term in the equation of motion for Q2. The first part of F is the usual 
stochastic force arising out of the fluctuations of the field in the vicinity of detector 2, while 
the second part is the retarded force due to detector 1. RSG correctly point out that these 
two forces are correlated. In the context of our formalism, these correlations are embodied 
in the correlation kernels 1/21 and z/12. 
Using the relation (4.2), we obtain 

F(r.) = ^ - ''-^P r <i--e-xi" r ds,,(s)e-'. (4.9) 

aT2 47r aT2 J-oo J-oo 

Consider the correlator of F with the correlator of 7J2 subtracted out. We have: 

< {F{T2),F{r^)} > < {^2{t2)M<)} >= (4.10) 

J —00 J —00 J-oo J —00 



TT^ dT2 dT2 
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The kernels i>2i and separate into advanced and retarded parts. For the advanced parts, 

6162 f°° dk 



^i\V2,s) = -^cosk{v2{r2) -vi{s)) 

6162 r°° dk , / / X 62 , , / . . . X 

= -rCos/c^^iTi^ =-^ii^fl,^ 4-11 

47r Jo k e\ 



and similarly 

DUs.T',) = ''-^i)l,{s,T'^). (4.12) 

ei 

The advanced parts of the correlations can therefore be constructed from the advanced part 
of the noise along the trajectory of detector 1. With this simplification, we obtain 

< {F(r2),F(r^)} > < {rj,{T,) Mr'.)} >= 

-T--rj[ dujujxl' dse |e «— z/h^tr, sj + e "td-^^hIs, r^jj 

TT aT-2 J -00 J -00 UTr QTr 

-i^ r du r da;We''"^«e^'^'^«xi'^xiV H ds H ds's-^'e-'""''' i>n{s, s')] 

TT J —00 J— 00 J— 00 J— 00 

+ {r.p.) (4.13) 
where (r.p.) denotes the retarded part: 



(r.p.) = (4.14) 

TT J-00 J-00 dT2 dT2 dT2 0X2 

At this point, we specialize to the case when detector 1 is uniformly accelerated. Then 
we have 

vi{n) = a-^e"-^'; ui{n) = -a^^e^^^^ (4.15) 

As shown in subsection 3.2, the noise i>ii is thermal and isotropic. The retarded time tr — 
ln(af2(T2)). We may substitute for i/n in equation (4.13) and carry out the integrations 
over s and s' to obtain 

< {F(r2),F(r^)} > < {v2{r2)M<)} >= (4.16) 

2hra2drndr'^ r^,^,_^,^nk 



, ,f r dkk coth - cos k{rn - r;,)(xi'^ + X?^* - 47i | x?^ H + (r.p-)- 
TT aT2 aT2 Jo a 

(1) (D* 

The first term in the above expression vanishes as a consequence of the identity Xk +Xk — 
471 I X^k^ mentioned earlier. The only contribution to the excitation of the probe is 
therefore from the retarded parts of the correlations z/12 and z/21. This asymmetry between 
retarded and advanced parts is really a consequence of the choice of retarded boundary 
conditions in the formulation of the problem (the states of detector and field are assumed to 
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be uncorrelated at past infinity) and the switching process at U2 — 0. The vanishing of the 
first term in the above expression is a gencraUzation of the cancellation obtained by RSG 
for a probe moving along an inertial trajectory. In order to study the retarded contribution 
in greater detail, it is desirable to simplify the correlations i>[2 and i>2i- The functions 
and Z21 take the form 

4:71 Jo k 

= ^ r^e-^(^-^+-(-)). (4.18) 
47r ^0 k 

Introducing the Fourier transforms 

e-^ ' ^ — dLue-'^^'r{ )(-)ve2^, k>0 

ZTTO 7-00 a a 

and substituting in (4.19), we obtain 

Z{.M) = da;e--ei^r(-^) T ^^(^)f e^^-(^^). (4.20) 

OTT O J—oo OL Jo K Qi 

The integral over k can be evaluated in terms of gamma functions. Explicitly making use of 
the fact that U2{t!^ > 0, and simplifying, 

^i2(n,r2) = ^y • (4.21) 

a 

Similarly, 

7r/ /^ eie2 [-^ duj COS ^{aT[ + lu aU2{T2)) 

a 

Therefore there is no contribution from the imaginary parts of and Z21 for U2{t!^ > 0, 
as obtained earlier (equation (4.5) leading to (4.6)). 

Differentiating the above expressions with respect to T2 and r2, and substituting in the 
expression for (r.p.), one obtains, after carrying out the integration over s, 

< {F{t2\F{t'2)} > < {V2{T2)Mr'2)} >^ M (4.23) 



2/i7i72 dkk (1) 2 / /N / ..ik_idu2dv2 ,2 / n / ,..ik_^du2dv; 

^ L ^i^^^ Mr2Mr2)) ^ + (« -2(r2).2(r2)) ^ 



2i 



2 L<,/2 dT2 dT2 



The coincidence limit of the above expression yields the fluctuations of the random force 
acting on the probe: defining 5F{t) = F{t) — ^?72(t), we obtain 

^„9, , 2h'yi'y2 du2 dv2 f°° dkk (d, n , / 1 
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The fluctuations are thus suppressed in the hmit of large U2V2 = — x^. For a probe 
trajectory without horizons, this is the limit in which the probe trajectory approaches future 
timelike infinity, which verifies that the effect of the accelerated oscillator on the field is 
ascribed to polarization rather than radiation (see also 0]). A radiation field is expected 
to persist at future infinity. 

Let us now turn to the question of the response of the probe. To obtain this, we will need 
to specify a particular form of trajectory for the probe as well. We will consider the simple 
inertial trajectory X2(t2) = 0, t2('i"2) = ^2, switched on at T2 = 0. Then equation (4.23) gives 

< {F(r2),F(r^)} > < Mr^),^^ >= (4.25) 

a 

Owing to the switching process at r2 = 0, the relation between Q2 and F (the Fourier 
transforms of Q2 and F) is a non-local one in frequency space, because of transient effects. 
However, if we restrict our attention to the late time behavior of detector 2, we obtain from 
equation (4.7) a local relation of the form 

Q2{uj) = ^F{u) (4.26) 

lUJ 

where 

F{u) = ^ r dre~'-F{T). (4.27) 
27r JO 

In the above expression, the lower limit of the r integration is zero, corresponding to the 
step function 9{t) multiplying F which enforces the switching condition. 
The correlator of Q2 is therefore given by 

< {Q2(^), g2(^')} >= < {F{u;), F{u')} > . (4.28) 

We have already obtained the difference of the correlator of F from its value in the absence 
of the accelerating detector, 1. Thus we have 

< {Q2{uj),Q2{uj')} >-< {Q^2\uj),Q^^\uj')} >= (4.29) 
L L dr'e-'--e-^-[< {F(r),F(r')} > --— < {r,2{r)Mr')} >] 



An^ujuj' Jo 

where the superscript (0) on Q2 refers to its value in the absence of the accelerating detector. 
Performing the integrations over r and r', we obtain 

< {Q2(o;),Q2(^0} > - < {Qf (^),Qf (^0} >= 

^7i72xi^^x2^ dkk (1) 2/^^x I i-i^-i / ^ , 



uu' J-oosinh^^'' 'a'' 



■nk 



{e^e{uo') + e-^e{-uj')). (4.30) 
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The step functions which distinguish positive and negative frequencies in the above expres- 
sion are an artefact of the switching process. 



5 Fluctuation-dissipation and correlation-propagation 
relations 

In this section we construct the fluctuation-dissipation relations for the detector system and 
extend this construction to obtain a new set of relations, which we call the correlation- 
propagation relations for trajectories without event horizons. These relations are a simple 
consequence of the analytic properties of the massless free field two-point function. We also 
discuss these relations in the context of the model of a uniformly accelerated detector and 
probe. 

Consider first the fluctuation dissipation relation for a quantum Brownian particle in a 
heat bath This can be expressed as a linear, non-local relation between the noise z/(s) 
and dissipation jl{s) kernels. Defining 7 by 

Ks) = J(^), (5.1) 
the finite temperature fluctuation dissipation relation is 

roo 

i>{s) = / ds'K{s - s')^{s') (5.2) 



where 

K{s - s') = coth(^) cos k{s - s') (5.3) 

JO TT 2 

is a universal kernel, independent of the spectral density of the bath. In particular, the kernel 
K is independent of the coupling constant e. Such a fluctuation-dissipation relation holds 
for the uniformly accelerated detector (with temperature given by the Unruh temperature) 
and the inertial detector (with zero temperature). It was derived in [^] in the context of 
a quantum Brownian model with bilinear coupling between bath (field) and particle (detec- 
tor). In such a model 7 is indeed the quantity which characterizes dissipation in the effective 
Langevin equation for the particle. 

In the context of the minimally coupled model, however, we find it suitable to define 7 

as 

7(3) = -^m (5.4) 

as 

as this is the quantity which directly appears in the dissipative term of the Langevin equation 
derived above. The fluctuation-dissipation relation then takes the form (5.3) with K defined 
as 

Kis-s') = r ^coth{(^) cos k(s-s'). (5.5) 
JO nk 2 
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An important aspect of either form of the fluctuation-dissipation relation is that the noise 
and dissipation kernels, and consequently are stationary, i.e. they are functions of s — s' 
alone. 

We wish to investigate whether a suitable generalization of the above relation holds for the 
full A'"-detector system. To this end, we assume that the detector trajectories are everywhere 
timelike and consider first only the kernels as they characterize noise and dissipation in 
the dynamics of the detectors i. Wc also assume that the detectors are switched on forever, 
thus excluding transient effects due to the switching process. Using advanced and retarded 
null coordinates introduced earlier, we define 



denoting the advanced and retarded parts of the kernel 7. 

The timelike property of the trajectories implies that | ^ |< 1. Together with the fact 
that ti{Ti) are increasing functions of Tj, this implies that ^ and ^ are necessarily positive. 
It also implies that the functions Mj(rj) and Vi{Ti) have unique inverses, if they exist. This 
can be proved by way of contradiction: assume that Ui{Ti) = Ui{Tl) for some 7^ r/. Then 
we have Xi{Ti) — Xi{Tl) = ti{Ti) — ti{T-), which means that the points and r/ have light- 
like separation. This contradicts the fact that the trajectory is everywhere timelike. The 
uniqueness of is shown in the same way. 

These two properties lead to the following simplification in the expression for 'ju: 

Un:ri)^p{n-Tl). (5.7) 

Thus we see that, for an arbitrary trajectory, the dissipation or radiation reaction kernel has 
the same form and is always local. This fact has been used in obtaining the dissipative term 
in the equations of motion for the accelerated detector and probe (4.1 and 4.7). 
The fluctuation-dissipation relation now follows in a straightforward manner: 

/oo 
dsKi{Ti,s)^u{s,Tl) (5.8) 
-00 

where 



■^/2.(r.,r;) 



dm 



(5.6) 



roo (^/j 

= / — r[cOS k{Vi{Ti) -Vi{s))+ COS k{Ui{Ti) -Ui{s))]. (5.9) 

JO Zttk 

Wc now ask whether a similar relation holds between the real and imaginary parts of Zij, 
i ^ j- This would not be a fluctuation dissipation relation in the usual sense, as the real 
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part of Zij describes correlations of the field between points on different trajectories rather 
than fluctuations, and its imaginary part describes the propagation of radiation between 
one detector and the other, rather than dissipation. We will call such relations "correlation- 
propagation" relations. 

If points on different trajectories have space like separations, the relevant (defined as 
~^^) vanish as a consequence of the vanishing of the commutator of a free field for 
points at spacelike separations. This is simply an expression of causality in the detector 
dynamics. However, the corresponding correlation Vij need not vanish, and hence there 
cannot be a general relation between these two kernels. Such a situation is realized most 
clearly, for example, in the case of two uniformly accelerating detectors, one in the right 
and the other in the left Rindler wedge. The trajectories, although individually timehke, 
are spacelike separated everywhere. The corresponding 712 and 721 will therefore vanish 
identically. However, z/12 and v^x will remain non-zero, reflecting the highly correlated nature 
of the Minkowski vacuum state. 

If, however, none of the detector trajectories possess past or future horizons (in Minkowski 
space this is true in particular for geodesic trajectories, but not only for geodesic trajectories), 
then each of them will lie completely within the causal future of the others. In that case, we 
can obtain correlation-propagation relations relating separately the advanced and retarded 
correlations to their "propagating" counterparts. These relations follow from the fluctua- 
tion dissipation relations along single trajectories derived above, essentially by a method of 
geometric construction : defining 7^"- = —-^ and similarly Yiji we have 

7S(^^'^;) = ""-^K^n) - vM))^^. (5.10) 

Since the trajectory i does not possess horizons, the null coordinates Ui and Vi range from 
—00 to 00. Thus the functions ViV^^ and UiU~^ are identity functions over the entire real 
line. Then we obtain, similar to equation (5.6), 

7--(r.,r;) = ^S{t, - v-\v,{t;))) 

= ^%ir^,v■'{vJ{r;^)) (5.11) 

and 

%(r.r;^ = ^%in,ur'{u,{T^)). (5.12) 
The correlations may be constructed from the noises Da in an identical manner: 

^tiiri,r;) = '-^ffcosk{v,{n)-v,{vi\j{r',))) 

= ^i>uir^,v-\{T;^)) (5.13) 

where we have inserted the identity function ViV^^ in the first step. Also, 

^(^,r;) = %(7,,^.-S(^))■ (5.14) 
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These two sets of constructions for the propagation and correlation kernels in terms of the 
dissipation and noise kernels enables us to write down the correlation-propagation relations 
simply by invoking the fluctuation-dissipation relations (5.8) as they separately apply to the 
advanced and retarded parts of the noise and dissipation along single trajectories: 

/oo 
dsKf'''{Ti, s)7^"/(s, rj), (5.15) 
-oo 

i^f and being defined earlier (5.9). Since the quantities 'jij are really just 5-functions 
and the quantities Kf'^ are proportional to , these relations can be equivalently viewed 
as constructions of the correlations % from the noises Da. 

The above relations hold for trajectories without event horizons. In the example of the 
uniformly accelerated detector and probe, the uniformly accelerated detector trajectory does 
possess event horizons. This manifests in the property that the range of Ui is restricted to 
(—00,0) and the range of Vi to (0, oo). The probe trajectory, on the other hand, will be 
chosen to be free of horizons. We will also now assume that the probe is switched on forever. 
Then we can construct the correlations i>2i and the quantities 721 from i>22 and 722 exactly 
as described above, and obtain the corresponding correlation-propagation relations: 

/•oo 



V. 
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/oo 
dsKr{T2,s)r2({s,T[). (5.16) 
-00 



This simply follows by invoking the fluctuation-dissipation relation along the probe trajec- 
tory, as described above. However, it is of greater interest to know whether such relations 
would follow from the fluctuation-dissipation relation along the uniformly accelerated tra- 
jectory. As explained, this will not be completely possible because the accelerated trajectory 
possesses horizons. This difficulty shows up when one tries to write down a relation of the 
form (5.16) for the quantities 1)12 and 712. To do this, we first express the functions Zij in a 
different form. This was done in Section 4 (see the steps leading from 3.54 to 3.58) for the 
restricted case U2{t2) > 0, V2{t2) > 0. If we remove this restriction, we find 

ZUriA) = ^ T— [coth(— )cos-(aTi-ln|a^2(T^)|)^(^2(T^)) (5.17) 
47r io a; a a 



^12' 



^1,^2 = 4^/ — coth( — ) cos - an + In aM2 ^2 ^ -M2 T-^ (5.18 
An Jo uj a a 

cos -(ari + In I aM2(T9) I) , , uj , , , ,^ ,^ , , 

+ sinhB^ '^''^' eiu2{T^))-isin-{an + \n \ au2{T^) \)0{-U2{t!,))]. 

a 

Z21 can be expressed in a similar way. From the above, we see that the advanced (retarded) 
correlation for V2 > {u2 < 0) has a thermal form, because these correlations can be 
constructed simply from the noise along the accelerated trajectory. We are therefore able to 
write down a correlation- propagation relation for this part of the correlations alone. This 
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takes the form 



/oo 
dsK,{n,s)^i2{s,T^) (5.19) 
-oo 



where 



luM) = (5.20) 



dTi 

6162 



[(5(ti - a-Mn | a^;2(r^) |)^(t'2(r^)) + 5{ti + a~Mn | aM2(r^) |)^(-M2(r^))] 



and 



ii'i(Ti,s)= / — coth — coskiTi — s). (5-21) 
70 tt/c a 

The single relation (5.19), as opposed to separate relations between the advanced and re- 
tarded parts, is a consequence of the fact that the thermal noise is isotropic and therefore 
contains equal contributions from advanced and retarded parts. 

In the context of the analysis of section 4, where the probe is switched on at U2 = 0, 
(5.19) becomes 

/oo 
dsK1{r,,s)rAs.<)- (5.22) 



Viewed as a construction of 1)12 from z/^^^, this relation lies at the heart of the RSG cancellation 
in (4.16). Viewed alternatively as an extension of the thermal fluctuation-dissipation relation 
on the uniformly accelerated trajectory, it thus places the role of thermal equilibrium in the 
RSG cancellation on firmer ground. 

We now turn to the part of the correlations which do not partake in the correlation- 
propagation relation above. These are the advanced (retarded) correlations for V2 < Q 
{u2 > 0), containing the sinh~^ factors, and are not expressible in terms of the noise along 
the accelerated trajectory. Rather, they represent true correlations across the future (past) 
horizon. If we specialize to the case 1*2 > as in Section 4, then these are exactly the 
correlations which contribute to the excitation of the probe in the guise of (r.p.), equation 
(4.23). The probe may therefore be said to be excited by free field correlations across the 
future horizon. 

If we specialize to the simple probe trajectory a;2(r2) = 0, ^2(^2) = T2, then we have 
'U'2{t2) = V2{t2) — T2 and the expressions (5.18) for Z12 acquire a symmetric form. In this 
special case, we can write down a correlation- propagation relation for the entire kernel z/12, 
by relating the advanced part of the correlations across the horizon to the retarded part of 
the propagation kernel, and vice-versa: we then have 

7i2(Ti,T^) = ^[S{n - a-Mn I ar^ |)^(t^) + Sin + a'Hn \ ar^ |)^(-t^)]. (5.23) 

Choosing 

,^ f°° dk , , tt/c cos A;(ri + r{), 
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we obtain _ 

/oo 
dsK[{n,s)^u{s,T^) (5.25) 
-oo 

as a correlation-propagation relation in this special case. The above relation cannot be geo- 
metrically constructed from the fluctuation-dissipation relation along the single accelerated 
trajectory. So far, we have not been able to show the existence of such relations in more 
general cases. The extra piece in the interpolating kernel K[ comes from correlations across 
the horizon, as explained earlier. 



6 Summary and discussion 

To summarize our work and findings, we have presented a general formalism to treat an 
arbitrary number of detectors modelled as oscillators in arbitrary kinematic states, and 
minimally coupled to a massless scalar field in 1 + 1 dimensions. In this approach, the 
scalar field has been integrated out and the detector dynamics is described by a reduced set 
of effective semiclassical stochastic equations. These equations nonetheless contain the full 
quantum dynamics of the field. Our treatment can be extended to massive fields and higher 
dimensions by making appropriate changes in the two-point functions Zij. 

We studied four examples, starting with a single inertial and uniformly accelerated detec- 
tor, mainly to illustrate the new description, and culminating in the treatment of a uniformly 
accelerated oscillator and a second oscillator which functions as a probe. We show that there 
exist fluctuation-dissipation relations relating the fluctuations of the stochastic forces on the 
detectors to the dissipative forces. We discover a related set of correlation-propagation re- 
lations between the correlations of stochastic forces on different detectors and the retarded 
and advanced parts of the radiation mediated by them. 

In the analysis of two inertial detectors, we find that the change in the state of the field 
due to the coupling with either detector modifies the impedance functions of both detectors, 
and hence their dissipative properties. Also, this coupling introduces a mutual impedance 
which describes the change in the response of one detector due to the fluctuations of the 
field in the vicinity of the other one. The field fluctuations (noise) in this case are relatively 
trivial, and non-trivial effects can be ascribed mainly to the impedance functions. 

In the case of the accelerated detector and the probe, on the other hand, the noise due 
to field fluctuations and the field correlations between the two trajectories play a dominant 
role. Since the probe cannot causally influence the accelerated detector, the dissipative fea- 
tures of this problem are relatively trivial. Here, we find that most of the terms contributing 
to the response of the probe cancel out, leaving behind a contribution that arises purely 
from field correlations across the horizon. This cancellation was earlier pointed out |TT], |rB[ 
to be a consequence of the identity Xfc + Xfc = 47 | Xfc P or variations thereof, which is a 
form of fluctuation-dissipation relation. Although we utilize this identity in our calculation, 
we observe, however, that this really follows from the dissipative properties of the acceler- 
ated detector and its free uncoupled dynamics. It therefore does not explicitly involve the 
fluctuations of the field. We point out that this cancellation can instead be understood to 
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follow because the correlations between the accelerated detector and probe trajectories can 
be expressed partly in terms of the noise or field fluctuations along the accelerated trajectory 
alone, and also because of the isotropy of this noise. The expression of correlation in terms of 
noise can be equivalently viewed as a consequence of the correlation-propagation relations we 
obtain in Section 4, which are appropriate extensions of a generalized fluctuation-dissipation 
relation directly relating fleld fluctuations to dissipative properties. 

A distinct feature of the influence functional formalism as used in this paper is the 
assumption of an uncorrelated fleld-oscillator initial state. As argued in the appendix, an 
uncorrelated initial state is more readily realizable in the derivative coupling model. However, 
since the minimal and derivative coupling models are dynamically equivalent, we expect our 
flnal results to be essentially unchanged, in particular the results of detector response in the 
various cases studied. The discussion of fluctuation-dissipation relations can be reformulated 
as well in a way suitable to the derivative coupling model. 

We would like to mention possible extensions of this work to other problems. In discus- 
sions of the quantum equivalence principle [^, [1^, one compares the response of a detector 
moving on a geodesic trajectory in Minkowski space, and coupled to a quantum fleld, to its 
response along a geodesic of a spacetime with a homogenous gravitational fleld. The idea is 
to derive a suitable transformation on the state of the quantum fleld which yields the same 
detector response in both cases. If one can flnd such a transformation, the equality of the 
detector response in both cases constitutes a test of the validity of the quantum equivalence 
principle for local physical processes. However, a homogenous gravitational fleld deflnes a 
global inertial frame, and so one is inclined to believe that the equivalence principle would 
hold for non-local processes as well, such as the effective dynamics of two spatially separated 
detectors coupled to the same quantum fleld. We plan to investigate this and related is- 
sues, especially the implications of our flndings on black hole backreaction and information 
problems in later works. 



A Appendix 

A.l Infrared Problems with the Minimal Coupling Model 

The infra-red effects of the minimal coupling model are not trivial. In fact, as we will show, 
the proper treatment of this simple model, including an ultra-violet but not an infra-red 
cut-off, leads us to identify a super- select ion rule which prefers a particular class of bases 
in the model's Hilbert space. Using the minimal coupling model and this preferred class of 
bases is equivalent to a derivative coupling model used by Unruh and Zurek [UZ] and a basis 
of direct products of unperturbed fleld and oscillator states. 
The MC Hamiltonian is 

Hmc = H^o + ^iP- + f ■ (A.l) 

It is straightforward to show that the expectation value of Hmc for soft photon states (i.e., 
low energy eigenstates of the free fleld Hamiltonian) has a contribution proportional to the 
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inverse of their unperturbed energy. This suggests that the true low energy states of this 
model must have strong correlations between the field and the oscillator. If we reject an 
infra-red cut-off as unphysical, then we must conclude that states that are direct products 
of field and oscillator states will actually have energies very much higher than that of the 
true ground state of the model. 

The poor behaviour of the basis of unperturbed {i.e. e ^ 0) energy eigenstates reflects 
the fact that, since there are field modes at all frequencies, expanding around e = correctly 
requires degenerate perturbation theory. Since we are trying to set up an open quantum 
system, however, our choice of basis is not merely a matter of convention. Different bases 
can imply different partitions of the complete Hilbert space into 'system' and 'environment' 
subspaces. In the influence functional formalism, one traces over the final states of the 
environment, and assumes an initial state which is often a direct product of system and 
environment. If we change what we mean by 'system' and 'environment', the final trace 
becomes a different operation, and the initial state becomes a different state. 

It is not the model itself, of course, but only the naive basis that is badly behaved: the full 
Hamiltonian is quadratic, and hence equivalent to a set of decoupled harmonic oscillators. To 
understand the problems with the basis of unperturbed energy eigenstates, and to identify 
a better basis, we should diagonalize the full Hamiltonian. 

The MC Hamiltonian may be diagonalized by defining new creation and annihilation 
operators. If the original field and conjugate momentum operators are and tt{x), the 
diagonalizing annihilation operators are 

ttk = I ^ ( r dx[C{k)Ek{x)+iiimkx\[\k\(t){x)+ni{x)]] +eC{k)[iQ - \k\P] , (A.2) 
^2h-n\k\ ^ 

where Ek{x) = (1 - ■^) coskx-^ sin A;|x| and C{k) = [(1 - ■^)^ + ^]~^, for a massless 
field. 

Note that the set of field-like operators {ak,0'k} diagonalizes the entire Hamiltonian. 
There is no normal mode of the coupled system which corresponds even weakly to the 
unperturbed oscillator, and all of the normal modes contain very non-local excitations of 
the field. If we wish to consider the oscillator as an open system coupled to the field as an 
unobserved environment, then, the basis of exact energy eigenstates of the combined system 
will not be particularly convenient. The fact that this model is easily solved exactly does 
not make the system-environment problem trivial. 

We can now, however, determine the effect on the true ground state of the unperturbed 
oscillator raising and lowering operators. We find that either of these operators maps the 
ground state onto a highly excited state, whose expected energy is infra-red divergent. Con- 
sequently, observations that are restricted to the oscillator sector alone (as it is defined from 
Equation ( [A.l|) may be said to require infinite amounts of energy. Considering the field that 



appears in Equation (|A.1|) to be an unoberved environment is therefore unphysical. 

We can, however, change our basis so that the Hamiltonian appears more benign in 
terms of the transformed operators. We may effect this transformation using the unitary 
operator U = exp —^eQ(j){0). This transformation mixes the field and oscillator sectors. 
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and so changes what is meant by an observation of the oscillator alone. We can check that 
the new oscillator raising and lowering operators, acting on the ground state, now produce 
states whose energy is ultraviolet divergent, instead of infra-red. A physically plausible UV 
cut-off then renders this energy finite, and it becomes reasonable to consider the new field 
sector as unobserved. Furthermore, a direct product of the unperturbed ground state of the 
transformed field and any finite-energy oscillator state now has finite energy, and so is not 
unreasonable as an initial state for the coupled system. 

An alternative way of expressing the advantage of the transformed variables is to say that 
in order for a degree of freedom in the theory to be observable in isolation, it must require 
finite energy to excite or de-excite it without affecting other degrees of freedom, and it must 
also be spatially local (except possibly at small scales) . Such an observed degree of freedom 
will be some linear combination of the exact normal modes. If we require the co-efficients 
in this linear combination to vanish at high energies and remain finite everwhere else, we 
can ensure finite energy observability. If we require that the co-efficients are constant at low 
energies, we also ensure local observability. 

We would then like to have a basis in which this observed hnear combination "looks like" 
a harmonic oscillator coupled to a scalar field. Given our conditions of local and finite energy 
observability, our original basis does not provide this feature, but our second basis does. 

The transformed Hamiltonian Hmc becomes, in the new operators, precisely the Hamilto- 
nian of the Unruh-Zurek model. We have therefore found that, even if we begin by analysing 
the minimally-coupled model, we may be compelled in the end to study the Unruh-Zurek 
model (with a UV cut-off) instead. More convenient for our subsequent calculations than 
the Hamiltonian for this model is its Lagrangian: 

Ldc = \ {mQ' - f^Q') + \ j dk - k^%^k - -^Qi'^k + K) ■ (A.3) 

Here is the time-dependent, spatial Fourier transform of the field 0. 

Note that, even in the more benign Unruh-Zurek basis, distinguishing the oscillator as a 

system observable independently of the field directly implies that there must be a UV cut-off. 
One often argues that a cut-off is appropriate because one is not interested in accurately 
describing physics at inaccessible energy scales; but in the case of the oscillator coupled to 
a field, there must really be a cut-off in the coupling in order for there to be any accessible 
energy scales! 

A. 2 Correspondence between MC and derivative coupling models 

In the minimal coupling model, the derivative of the oscillator coordinate couples to the field, 
whereas in the derivative coupling model the oscillator coordinate couples to the derivative of 
the field. These two models thus differ by a total derivative term in the Lagrangian. In par- 
ticular, they have the same Heisenberg operator dynamics. The above subsection describes 
the issue of the initial state, and argues that an uncorrelated initial state is physically more 
realistic in the derivative coupling model. However, since the two models have the same 
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dynamics, this should translate to a simple prescription for switching from one model to the 
other in the context of the influence functional treatment. 

In the previous sections, we have derived all results from the minimal coupling model. 
One can obtain corresponding quantities in the influence functional of the derivative coupling 
model via the prescription below: 

dQi ^ 



dQ[ 



dn 

d d 



dTi dr'j 

The stochastic effective action in the derivative coupling model is then given by 



/I. [1 ) 
^ dn[Q-Qt - ^-QTQt - Q7 Siin)rjiin) 

-2Q-s,{n) f: ^'^"^^ dT'^Qpj{T'^)fl,j{n, rj)]. (A.5) 

-•1 J —oo 



Note that the quantities flij in the above equation refer to the newly defined quantities in the 
derivative coupling model. They are obtained by differentiating the corresponding quantities 
in the MC model twice. 

The Langevin equations are: 

-f^ - 2 ^ / ' dT'.Sj{T'.)si{Ti)jlij(Ti, r'AQj + QfQi = Si{Ti)r]i{Ti). (A. 6) 

dTi j=i''-°° 

The noise kernel, as the correlator of rji and rjj, is also obtained by the corresponding noise 
kernel in the MC model by differentiating twice, according to the correspondence established 
above. 

The infrared divergent energy of the initially uncorrclatcd state docs have an effect: the 
propagation kernel, in the MC model, contains an initial '"shock wave" term, as well as the 
expected dissipation and propagation terms, and this term is not present in the DC model. 
Since this shock wave is a transient, it has no significance in our late-time analysis. 
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